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We introduce a new approach to the periodic Anderson model (PAM) that allows a de- 
tailed investigation of the magnetic properties in the Kondo as well as the intermediate 
valence regime. Our method is based on an exact mapping of the PAM onto an effective 
medium strong-coupling Hubbard model. For the latter, the so-called spectral density ap- 
proach (SDA) is rather well motivated since it is based on exact results in the strong coupling 
limit. Besides the T = phase diagram, magnetization curves and Curie temperatures are 
presented and discussed with help of temperature-dependent quasiparticle densities of state. 
In the intermediate valence regime, the hybridization gap plays a major role in determin- 
ing the magnetic behaviour. Furthermore, our results indicate that ferromagnetism in this 
parameter regime is not induced by an effective spin-spin interaction between the localized 
levels mediated by conduction electrons as it is the case in the Kondo regime. The magnetic 
ordering is rather a single band effect within an effective /-band. 

PACS: 71.10.Fd, 71.28-hd, 75.30.Md 
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1 Introduction 

The lanthanides and actinides, and their compounds, 
show a great variety of interesting physical properties. 
Some of these have to be ascribed to the most unex- 
pected and least understood phenomena in condensed 
matter physics. Probably, the most diverse physical 
characteristics are found in the heavy-fermion (HF) and 
intermediate valence (IV) materials. A comprehensive 
review is given in ref . |Q| . 

Prototypes are the Gerium (Ge) and Uranium (U) in- 
termetallics. In these materials, incompletely filled inner 
4/-shells (in cerium, 5/ in uranium) are responsible for 
the unusual physical properties. Worth to mention are, 
of course, the name-giving heavy fcrmion phenomenum 
which is characterized by an enormous effective mass of 
crystal electrons and the anomalous superconductivity 
which is called anomalous due to its coexistence with 
magnetic ordering. The magnetic phase diagrams of 
these materials are quite extraordinary in their large 
variety of different phases, ranging from simple para- 
magnetic states to various kinds of magnetic ordering. 

Ferromagnetism is found in several Kondo lattice and 
HF materials. For instance, in CeSix phase transitions 
as functions of the Si concentration were observed 1^ , or 
in the heavy fermion compound CePd2Ga^, ferromag- 
netic ordering breaks down as function of external pres- 
sure For CeSix, the magnetic phase diagram with 
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respect to the Si concentration x seems quite clear, but 
some effects accompaying the magnetic ordering, like 
the irregular resistivity behaviour, are not yet under- 
stood. For the second material, CePd2Gaj,, not even 
the magnetic phase diagram is settled yet. The pressure- 
induced suppression of ferromagnetism could be followed 
by some other magnetically ordered phase. Besides of 
the above mentioned, a lot more materials showing fer- 
romagnetic ordering were found. 

Most attempts to understand HF- or IV-materials 
theoretically are based on the Anderson model Q . This 
model describes a system of uncorrelated conduction 
electrons which hybridizes with either a single localized 
electronic level (single impurity Anderson model, SIAM) 
or a lattice of localized levels (periodic Anderson model, 
PAM). 

Many approximation schemes for these models have 
been proposed. At least for the SIAM, there is enough 
confidence in some of these approximations that most 
of the results are widely accepted. The so-called non- 
crossing approximation (NGA)Q has proven to be in 
many aspects equivalent to quasi-exact QMG calcula- 
tions . In addition, by means of Bethe ansatz Q and 
renormalization group theories, many exact results could 
be obtained. 

For the PAM, the situation is not yet comparable. 
Even though many approximation schemes exist, begin- 
ning with the mean- field-approximation second or- 
der perturbation theory in J7 ^ or in the hybridization 
|Tc| , im , a modification of the non crossing approxima- 
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tion (LNCA) In2|, |3|, and various dynamic mean-field 
theories |15|, none of these are completely satisfac- 
tory in all aspects of the PAM. Some of the approxima- 
tions, namely the LNCA |l^ are restricted to the 
Kondo regime of the Anderson model. In this param- 
eter regime, denoted by a half-filled low-lying localized 
state, it is possible to map the PAM onto the s-f model 
in which the localized states are reduced to their spin 



degrees of freedom |16 



So, despite of all efforts done in this field, no fully sat- 
isfactory theoretical understanding of the 4/-materials is 
developed yet. Even on the experimental side, not every- 
thing has been clarifed. For instance, it is currently in- 
tensely discussed, whether the Kondo resonance should 
be seen in photoemission studies on periodic crystals 
||l7| , ^ But these disputs will probably be settled 
in the near future. 

For the theoretical side, it still seems neccessary to de- 
velope better theories which help to understand at least 
partial aspects of the PAM. 

In this paper, we will introduce a novel approx- 
imation scheme for the periodic Anderson model us- 
ing an exact mapping of the PAM onto an effective 
Hubbard model. The model parameters within this 
effective model will be such that results obtained in 
a strong-coupling perturbational theory are significant. 
These motivate the spectral density approach (SDA) 
1^, H to solve the effective Hubbard model. This 
approximation scheme has proven, at least for the Hub- 
bard model, to be trustworthy concerning the magnetic 
properties of the system. 

This approximation scheme will allow to solve the PAM 
for a wide range of system parameters, the crossover 
from the Kondo limit to the intermediate valence region 
will be subject of our investigation. Although due to 
numerical simplicity, all our results were obtained un- 
der the assumption of a fc-independent selfenergy, this 
method is in principle not restricted to the local approx- 
imation as dynamical mean-field approaches are. Fur- 
thermore, even if based on exact results obtained in the 
limit U — > oo, it is not of perturbative character, there- 
for not restricted to this limit. The major drawback is 
the neglection of quasiparticlc damping, and therewith 
the impossibility to show fermi liquid behaviour. 

We will develope our approximation scheme in the 
next section of the paper. In section three, the results 
obtained with this method concerning ferromagnetism 
in the periodic Anderson model are presented. 

2 Theory 



2.1 Model Hamiltonian and its many- 
body problem 

Starting point is the periodic Anderson Hamiltonian 



i,(T i.{7 

Si,cr ifi.a) ^ud s| ^ ificr) '"^^^ ^^'^ annihilation and cre- 
ation operators for an electron in a non-degenerate 
conduction band state (locahzed /-state), and nlj = 
fiafi,iy is occupation number operator for the /- 

states. The hopping integral J = ^Sge-*^(^'-'^^'e(fc) 
describes the propagation of free, i.e. unhybridized con- 
duction electrons from site j to site i, e/ is the position 
of the free /-level relative to the center of mass of the 
conduction band density of states. The hybridization V 
is taken as a real constant, and finally U is the Coulomb 
repulsion between two /-electrons on the same lattice 
site. By use of the commutators 



(2) 



[/,,,, H]_ - £//,, -f Vs,^, + UU^n^l^ (3) 
we find for the single-s,/ electron Green's functions 
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the following equations of motion: 



(6) 



{E-ei)iU,,-!l^))E= (7) 
M,,, -I- fl^E + UU,,4fl^; f]jE 

After Fourier transformation these equations read: 

{E - <k))G^i'jE) = h + y((/,, ^; s\ j)E (8) 
{E-ef)d^l{E)^ (9) 
^+y^%.^fljE + UD^,^^{E) 
The "higher" Green's function D^^{E), 

^kJE) = ^ E«/pV./?-4+fe-,>; 4 (10) 



2 



prevents a direct solution of the equations of motion. 
However, the "mixed" function in (||) and (^), respec- 
tively, can easily be eliminated using the respective 
equation of motion: 



(11) 



So we are left with the following set of equations of mo- 
tion: 



h+ '-^Gi^^ (E) 



V 

E-e{k) 



(12) 



E-e{k)' ' (13) 

^+UD^^JE) 



which tells us that the determination of the /-Green's 
function (jj) solves the problem. 

The introduction of a selfenergy Yjj: ^{E) by 



^kAE)GfiE)^UDj:jE) 



(14) 



allows a formal solution of the equations of motion ( p^ ) 
and (PI): 



I E - (E) - ef 
G (E) — h 

{E-J:f:jE)-ef){E-e{k))-V^ 

G^iE) = 



J E-e{k) > 



(15) 
(16) 



From these we derive the quasiparticle densities of states 
(QDOS): 

P^AE)^^Y.{-\^GfJE^^O^)) (17) 

k 

PfAE) = ^Y.{-\^Gi^y + ^^^)) (18) 



The full problem is then obviously solved as soon as 
we have found the selfenergy Eg^(£'), defined by (p^. 
Our approximation scheme is developed in the next two 
sections. 

2.2 Effective medium approach 

The basic idea of the approximation is a mapping of 
the periodic Anderson model onto an effective medium 
Hubbard model. This effective Hamiltonian will have 
an energy-dependent one-particle term, containing the 
influence of the hybridization of the localized states with 
the conduction band and a Hubbard-like Coulomb term 
for the respective quasiparticles and has the following 
structure: 



^(^)-E^^('?)«L«£,. + J^E-i: 



(a)^(a) 



with 



Afe(^) 



rj — e(k) 
1 



N 



(22) 

(23) 
(24) 

being the effective band dispersion and nl^^ — al^Ui^^. 
•q has to be understood as a parameter with the dimen- 
sion of energy. Throughout the paper, we will denote 
quantities which are defined in the effective medium, 
with a tilde (~). All these quantities depend on the 
parameter rj which will often be omitted for readabil- 
ity. The Green's function for the system described by 



ik(Ri-Rj) 



(|22|), g''^ (E) = {{at ; al )) e can be determined by the 



equation of motion: 



{E - Aj:irj)) &^>^{E) UOf^iE) (25) 



The higher Green's function (E) is defined analo- 

gously to ( [ic| ) and can be expressed in terms of a prop- 
erly defined selfenergy: 



In the case of a local selfenergy, 

Y.~^JE)^Y.^{E) 



(19) 



simple manipulations yield the following structures of 
the densities of states: 



PsAE) =Po ( E 
PlAE) = 



E-^,iE)~ef 

y2 



{E~i:,iE)~ef] 



X pqIE 



E-j:4E)-ej 



(20) 
(21) 



with po{E) — -i-SjrM (^E - e(fc)^ being the free s- 
electron Bloch density of states. 



eI") (E)G'y (E) = UD'y {E) 



■.in), 



(26) 



By comparision of the two equations of motion (|T^) and 
(25), it follows that 



r]=E 



and equivalently for the respective selfenergies 



f!^jE) 



r]=E 



^kJE) 



(27) 



(28) 



Thus, by solving the effective problem defined via (|22 
for all values of rj, one can obtain the /-Green's function 
of the periodic Anderson model. The original problem 
has been mapped onto an "effective" Hubbard problem. 
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Since the Hubbard model, at least in more than one 
dimension, is not yet solved exactly, the advantage of 
the mapping is not seen immediately. But a closer anal- 
ysis of the effective medium will show that the mapping 
allows for a rather well motivated approximation. 

The "free" band dispersion Aj:{ri) changes with 77 be- 
ing, however, for reasonable 1^-values always very nar- 
row. The respective "free" density of states 



k 



E-ef 



(29) 



is plotted in figure |l| for some typical examples of rj. 
The effective dispersion Ag(77) diverges for e(fc) = 77, 




_ Ti=0.7 
Ti=0.25 



Figure 1: Effective density of states for different values 
of 77. The free Bloch density of states is taken to be 
semielliptic of unit width and center of gravity at i? = 0. 
The hybridization is F = 0.2 and the /-level ef — —0.02 



i.e. when rj falls into the band region (mine(A:) < rj < 
maxe(fc)). The resulting density of states extends from 
—00 to 00. But since p^^\e) ^ the respective 
weight vanishes. Thus, the effective width W is very 
small scaling with the hybridization V which must be 
considered for all realistic situations as a small param- 
eter, at least very much smaller than the Coulomb in- 
teraction U . The "Hubbard problem" which is given by 
the effective Hamiltonian ( p^ ) is therefore to be ascribed 
to the strong couphng regime: 1. 

At this point, the advantage of the mapping onto 
the effective model is clear. We have replaced the peri- 
odic Anderson model, which is, though a minimal, still a 
two-band model with an effective single band Hubbard 
model, for which some standard approximation methods 
exist. By construction of the effective medium, the sin- 
gle band model turns out to belong to the strong cou- 
pling regime, thus results obtained in a perturba- 
tional theory ]2^ can help constructing an appropriate 
approximation scheme 



2.3 Spectral density approach 

We use a selfconsistent spectral density approach 
(SDA)pO|, pl| , p2| to find an approximate solution of the 
"effective" Hubbard problem defined by the Hamilto- 
nian H{vi) in (p^). The method is based on a physically 
motivated ansatz for the single-electron spectral density. 
Its main advantages are the physically simple concept 
and the non-perturbative character being not restricted 
to Fermi-systems but also working for Bose- and even 
classical systems ^ . Recent applications of the 

SDA concern the attractive {U < 0) Hubbard model 
p7| ] , the t — J model , the magnetism and electronic 
structure of systems of reduced dimensions as thin films 
and surfaces |3^, ^ . It is also used for the investiga- 
tion of high- Tc superconductivity 1 32 , 3^ . In this paper 



we apply the SDA to the "effective medium" -Hubbard 
model (p^). 

The single-electron spectral density is defined by 

s^,,AE)^^ I d£;e-*^*(K.W;«L(o)]+) (30) 



where [...;...]+ denotes the anticommutator and (...) 
the thermodynamic average. The construction opera- 
tors are taken to be in the Heisenberg time-dependent 
picture. 

In an exact spectral-moment analysis in the limit 
J7 — > 00, Harris and Lange have shown that the spectral 
density essentially consists of a two-peak structure [pSl . 



The effective Hubbard model (22) must be considered 
in the strong-coupling regime. Therefore it seems ap- 
propriate to make the following ansatz for the spectral 
density: 



k 



5{E 



(31) 



The still unknown parameters Ej^cy and 5^,0- , the quasi- 
particle energy and spectral weight, can be calculated 
by use of the first four moments of the spectral density: 



+00 



k,(T 



dEE"S-,jE)^ 



(32) 



-fold commutator 



This procedure is identical to the one performed in |22 
for the real Hubbard problem. An explicit description 
of the calculation is presented there. 

Despite its obvious restrictions, e.g. the complete 
neglection of quasiparticle damping, the two-pole ap- 
proximation together with the moment method to cal- 
culate the free parameters is able to describe the mag- 
netic properties of the Hubbard model surprisingly well 
|23,|34| 
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As a result one obtains a selfenergy of the following 
structure: 



E - - Ft 

° k, — <7 



U{1 (n^l)) 



(33) 



The decisive terms are and Fj^ _^ which distinguish 
this selfenergy from the Hubbard-I solution |3^. These 
terms, mainly consisting of higher correlation functions, 
may provoke a spin-dependent shift and/or deformation 
of the bands and may therefore be responsible for the 
existence of spontaneous magnetism ^ . 

The fc-dependent term Fj: _^ seems to be of mi- 
nor importance for the magnetic behaviour. Since 
SgFg _^ = 0, it does not change the center of grav- 
ity of the density of states being mainly responsible for 
a deformation and narrowing of the bands. We have 
neglected this term in the following calculations. A de- 
tailed inspection of the influence of this term, as well 
as quasiparticle damping will be the subject of future 
investigations. 

The term B^^ has the following structure: 



= An 



1 



)) 



■ct"'J,-o' 



i2n^ - 1)) 



+ CXD 



(34) 



U 



Aq denotes the center of gravity of the effective electron 
hopping Aij (see equation (p^). In the case of a k- 
independent, real selfenergy, as we have obtained with 
our approximation, the spectral density can be expressed 
as: 



Sj: jE) = hs[E- A,-(r?) + /i - ^(J'He) 
Together with the spectral theorem: 



(35) 



(36) 



the equations (|34D, (|33|) and (|35|) form a closed set of 
equations which has to be solved selfconsistently. 

After retransforming from the effective to the real 
medium according to equations ( p7| ) and (|2^) respec- 
tively, one finally obtains a solution for the periodic An- 
derson model. 



3 Results 

In the following sections, the results of the above devel- 
oped theory will be presented. The system is determined 
by the following parameters: The free conduction band 
is described by a semi-elliptic density of states of the 
width W = 1 and center of gravity at i? = 0, thus defin- 
ing the energy scale. The /-level is characterized by the 
parameters ef, its relative position with respect to the 
center of mass of the conduction band, and the on-site 
Coulomb interaction U. The hybridization V is taken 
as a relatively small parameter, in all cases examined 
below, we chose V = 0.2. Of crucial importance is the 
total number of electrons per lattice site, ntot = Uf + Us- 
Calculations will be made for the ground state {T — 0) 
and for finite temperatures. 

3.1 The paramagnetic system 

In figure ^ the paramagnetic quasiparticle density of 
states is plotted for different values of e/. The inter- 
action constant is chosen to be C/ = 1.5, the total occu- 
pation number is ntot = 1-5. The arrows point to the 
position of fi, the chemical potential. 




1.0 0.0 1.0 2.0 3.0 



Figure 2: s- (dashed line) and /(solid line)-density of 
states. The parameters are U = 1.5, V — 0.2, T = 
and ntot — 1-5. The chemical potential /i is indicated by 
the arrows. 



It can be clearly seen that the density of states is di- 
vided into three structures. With the /-level positioned 
within the free conduction band, the lower two struc- 
tures exhibit band-like behavior, while the third struc- 
ture, the upper "Hubbard-band" of the /-level, retains 
its atomic shape. 
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The gap which separates the lower two peaks, is in- 
duced by the hybridization. Its position is related to the 
position of the free /-level, and its size scales with the 
hybridization strength V. 

Contrary to the upper "Hubbard-band" which is al- 
most of pure /-character, these lower structures show 
strong mixing of s- and /-states. For approaching 
one of the band-edges of the free s-band, a strong peak 
of /-states develops. 

In the example shown in figure ^ there is always a 
peak close to the Fermi energy. This coincidence origi- 
nates from the chosen system parameters, more precisely 
from the number of electrons ritot- 

Due to the simple concept of the two-pole approxima- 
tion used for solving the effective Hubbard model, we do 
not see a Kondo resonance. A Kondo resonance would 
strongly influence the low energy properties of the sys- 
tem. But as already discussed in the theory section, for 
describing the magnetic properties, the correct high en- 
ergy behaviour of the selfenergy is of importance ^T\ . 
Therefore, we do believe in our results concerning mag- 
netism at least for situations and temperatures where a 
Kondo effect is not important. 

3.2 The magnetic ground state proper- 
ties 

In figure ^, we have drawn a magnetic phase diagram of 
the periodic Anderson model derived within the above 
developed theory. This phase diagram consists of an 
area where ferromagnetism is stable whereas in the re- 
maining parameter space (ritot and e/ being the vari- 
ables) the ground state is paramagnetic. Antiferromag- 
netic ordering was not considered. 



2.0 




0.0 I ■ ' ■ ' ■ ' ■ ' ■ 1 

-0.8 -0.6 -0.4 -0.2 0.0 0.2 

e, 

Figure 3: Magnetic phase diagram for the periodic An- 
derson model at T = with [/ = 10 and ^ = 0.2. The 
shaded region is ferromagnetic, the surrounding param- 
agnetic. 

For a fixed total number of electrons ritot ferromag- 



netism can only exist if the /-level lies in a certain 
range somewhere around the lower half of the conduc- 
tion band. For energetically lower positions of the /- 
level, the ferromagnetic solution of the set of equations 
still exists. But since the paramagnetic solution has a 
lower energy, the ground state is paramagnetic. With 
the /-level approaching higher energies relative to the 
conduction band, the number of electrons in /-states, 
(n/) reduces, and the /-density of states broadens. It is 
commonly accepted that both effects reduce the stability 
of ferromagnetism in the system. So it is not surprising 
to see the system become paramagnetic with higher val- 
ues of e/. 
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Figure 4: /- and s-magnetization and the occupation 
number Uf oi the /-level as functions of e/ for different 
total occupation numbers ritot- {U = 10, V ~ 0.2, T = 
0) 

In the following, we want to focus on the intermedi- 
ate valence regime (IV) of the periodic Anderson model. 
Intermediate valence is characterized by non-integer val- 
ues of n/, the /-level occupation number. Per defini- 
tionem, an "atomic-like" level can contain either none, 
one or two electrons. Due to the rather high value of 
U, the local part of the Coulomb repulsion, double oc- 
cupancy is suppressed. So a non-integer value of ritot 
describes mainly fluctuations between the non-occupied 
and the single-occupied state. The condition for IV is 
found for e / lying within the conduction band, more spe- 
cific, the closer e/ and n are, the higher are the charge 
fluctuations within the /-states, the average rif occupa- 
tion number is lower. 

For describing rare earth materials, a total occupa- 
tion number ritot lower than ritot = 1 is not of interest. 
So we will confine our investigations to the parameter 
range which builds up the upper right hand side of the 
phase diagram in figure ||. In this region, there is a "bay" 
in the ferromagnetic isle. For a total number of electrons 
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around 1.75, ferromagnetism first disappears with rising 
/-level, but reappears for little higher /-levels, until the 
system finally becomes paramagnetic again. 

This is shown in more detail in figure ^, presenting 
the magnetization as a function of e/ for three different 
electron densities Ufot = {1-5, 1.6, 1.7}. For Utat = 1.7 
the two separate regions of stable ferromagnetism can be 
seen, but already for the lower occupation numbers, two 
distinct regions of ferromagnetism can be distinguished 
from each other. 

The first of these, appearing for low values of e/, sets 
in with a first order phase transition at e'^^^^'^ and is 
stable up to e^'J.^j"'^. The existence of a first order phase 

transition seems unusual. But since e^™^'^ is a quan- 
tum critical point, not a thermodynamic critical point, 
first order phase transitions are not forbidden. At least 
for transitions from antiferromagnetic to paramagnetic 
states, first order phase transitions were experimentally 
observed Q. 

This region, usually called Kondo regime, is charac- 
terized by almost half filled /-levels forming relatively 
large local moments. According to the Schrieffer- Wolff- 
transformation fl^ ], one expects to find antiferromag- 
netic coupling between the /-moments and the conduc- 
tion band moments. This behaviour is clearly seen in 
figure ^. But the uncorrelated conduction band shows 
only a weak magnetic polarization. Thus, the indirect 
interaction between the local moments via the conduc- 
tion band is responsible for the magnetic ordering. This 
region of ferromagnetism could be called "local moment 
magnetism" (LMM) , since the electrons forming the rel- 
ativly large local moments are energetically far away 
from the Fermi energy. Thus, they do not participate 
in transport phenomena. 

The second region of ferromagnetic ordering belongs 
to the intermediate valence regime of the PAM. This is 
specified by /-level occupation numbers Uf away from 
half filling, and a typical broadening of the /-levels to 
bands. The average magnetic moment at each /-site is 
strongly reduced. Furthermore, the coupling between 
the /-levels and the conduction band is not always an- 
tiferromagnetic, but, depending on the system param- 
eter, can be ferromagnetic. So, intermediate valence 
magnetism (IVM) is in many aspects different from the 
LMM. 

The crossover from the LMM into the IVM regime 
depends on the electron density. As already men- 
tioned, for higher ritot, ferromagnetic order disappears 

, upper, 1 J , lower, 2 • j j i 

at e^^^j^ to reappear at ^^.j^. m a second order phase 
transition. For ntot = 1-6, the two critical points coin- 

• j . upper. 1 lower, 2 ,1 - j- i_ j_ 

cide. I.e. = ^/crit ' ^WLS an immediate transi- 

tion between the two regions of ferromagnetism is seen. 
For even lower occupation numbers, e.g. ntot = 1-5 there 
is a crossover region with linear behaviour of the mag- 
netization. 

The reason for these different kinds of crossover be- 
haviour can be found in the quasi particle densities of 
state (QDOS), as shown in figures | to |. These QDOS 
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Figure 5: /- (s-) quasiparticle density of states for dif- 
ferent values of e/ in the upper (lower) row (full line: 
spin-|, dotted line: spin-|). The arrow denotes the po- 
sition of the chemical potential. The total number of 
electrons is ntot = 1-5, the other system parameters are 
V = 0.2, T = 0, and J7 = 10. 
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Figure 6: Same as figure |^, but ntot — 1-6. 
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Figure 7: Same as figured, but ritot = 1-7. 



correspond to the three cases presented in figure ^ with 
ntot = {1-5, 1.6, 1.7}, respectively. In all three cases, 
the crossover region coincides with the region where the 
chemical potential fi crosses the hybridization gap. In 
case of fcrromagnetism, the hybridization gap of the ma- 
jority spin direction band is decisive. 

For ntot = 1.5 (figure ||), the chemical potential lies 
in the flat region of the minority band, thus a small 
change of e/ will, for constant ntot, have only minor ef- 
fects on the quasiparticle densities of states. 
At ritot = 1.6 (figure ||), the different behaviour is pro- 
voked by the fact that the spin-j band has a rather sharp 
peak around /i, so already small parameter changes can 
have a strong impact on the system. 
Finally, for ritot = 1.7 (figure fcrromagnetism breaks 
down. 

It is worth to mention the possibility of a metal- 
insulator transition due to the hybridization gap. For 
ntot = 1.7, we see such a transition for e/ « 0.1 as 
the system becomes paramagnetic. On the contrary, for 
ntot = 1.5 or utot = 1.6, the system stays ferromagnetic 
during the transition, and only spin-J,-electrons can con- 
tribute to conductivity. But no metal-insulator transi- 
tion is seen in this case. 

Another fairly striking result of our calculations is 
the possibility of ferromagnetic coupling between the 
conduction band and the /-states. 

The Schrieffer- Wolff transformation maps the PAM 
onto the Kondo model jl6| . This model considers only a 
spin-exchange interaction between the conduction band 
and the system of /-levels. The effective coupling be- 
tween the different electron systems is antiferromag- 
netic. So, one would expect to find the conduction band 
spin-polarized anti-parallel to the /-states. 

As already mentioned above, we do find parameter 
sets, where the conduction band and the localized lev- 
els are ferromagnetically coupled. This can be seen in 



figure ^, and in more detail in figure ||. There, the /- 
and s-magnetization as well as the /-occupation num- 
ber nf is plotted as a function of ntot for three different 
values of e/ all very close to the center of gravity of the 
free conduction band, thus belonging to the intermedi- 
ate valence regime. Dependent on ntot the conduction 
band is polarized either in the same spin direction as the 
/-system (ferromagnetic coupling), or in opposite direc- 
tion (antiferromagnetic coupling). For higher nj, the 
tendency is clearly towards antiferromagnetic coupling. 



0.70 



0.30 " 




1.30 1.40 1.50 1.60 1.70 1.80 

n,„, 

Figure 8: s- (/-) magnetization nig (nif) and the /- 
occupation n/ as a function of the total occupation num- 
ber ntot- In all four cases the free /-level e/ is close to the 
center of mass of the free conduction band. The other 
system parameters are U = 10, V = 0.2 and T = 0. 



If one keeps in mind that the Schrieffer- Wolff trans- 
formation is only valid in the Kondo limit of the PAM, 
i.e. for half filled /-levels, no discrepancy is found. 

The mechanism which controls the direction of spin- 
polarization of the conduction band can be seen in the 
densities of state (figures ^ to |ll|). While for pj^^iE) the 
weight in each of the quasiparticle bands is distributed 
very asymmetric, ps^a[E) is always rather flat. But due 
to the hybridization, the partial s-electron QDOS ex- 
pierences the same spin-dependent bandshift as the re- 
spective /-QDOS. So, especially, when the lower band 
is not yet filled (i.e. low ntot), the spin-t s-band out- 
weighs the respective spin-| band. Beginning with /i 
crossing the upper band edge of the lower band, the 
spin- J, s-band gets filled, thus the QDOS induced ten- 
dency towards ferromagnetic coupling gets reduced and 
the antiferromagnetic exchange dominates. 

The antiferromagnetic exchange expresses itself in 
Ps,a{E) by the fact that for spin t the lower of the hy- 
bridization bands has less, for spin J, more weight than 
the upper hybridization band. 

The change of sign of the conduction band polariza- 
tion provokes another interesting question. There is a 
dedicated set of parameters for which the s-band is not 
polarized at all, but the /-moments order ferromagnet- 
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Figure 9: /- (s-) quasiparticle density of states for dif- 
ferent values of ntot in the upper (fower) row (full line: 
spin-f, dotted line: spin-j). The arrow denotes the po- 
sition of /i, the chemical potential. The free /-level is 
located at e/ = —0.02, the other system parameters are 
V = 0.2, T = 0, and t/ = 10. 
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Figure 10: Same as figure ^, but with e/ = 0.00. 
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Figure 11: Same as figure O, but with e/ = 0.02. 



ically. This suggests that the spin polarization of the 
s-system is not responsible for the magnetic ordering of 
the /-levels. 

In the LMM regime, the situation is clear. A strong 
antiferromagnetic spin exchange between localized levels 
and conduction band states leads to an indirect coupling 
of the /-moments, thus to ferromagntic ordering. 

On the contrary, in the IVM regime, ferromagnetism 
in the /-system is obviously independent of a magnetic 
ordering in the conduction band. Otherwise one would 
expect to see some effects on the /-magnetization at the 
critical occupation number ntot where = 0. In figure 
^ nothing like this can be observed. So in our opinion, 
spin exchange processes between the conduction band 
and the localized levels cannot be the major effect for 
producing the ferromagnetic ordering in the /-system. 
Of course, spin fluctuations are still present, thus induc- 
ing the magnetic ordering in the conduction band. But 
there must be a different mechanism leading to the fer- 
romagnetic ordering of the /-moments. We propose a 
single band effect within an effective /-band. Due to 
charge fluctuations between /- and s-system, an effec- 
tive /-/-hopping is possible. This leads to the formation 
of a narrow band analogously to a single-band Hubbard 
model. In the latter, ferromagnetism is, at least un- 
der certain conditions, expected to exist [|9[ p2[ , 
favored by a rather small band which should show a 
divergence at a band edge, and an occupation number 
not too far away from half filling. By examining the /- 
electron QDOS (see for example figure ^, one sees that 
these conditions are fuUfiUed in the effective /-band. 

3.3 Finite Temperatures 

When examining magnetism, the temperature depen- 
dence of several quantities is of major interest. In this 
section, the magnetization curves as well as the Curie- 
temperatures Tc for different system parameters are pre- 
sented. 

When the bandwidth of the free conduction band 
is taken to be W = 1, the temperatures are given in 
Kelvin. The figures |l^ and |l^ show the temperature de- 
pendence of the s- and /-magnetization for two electron 
densities and various /-level positions e/. The T — 0- 
properties for these parameters were presented in figure 

For Ufot = 1.5, as presented in figure 03, the magne- 
tization vanishes continuously at Tc, whereas for ntot — 
1.7 (see fig. |l3|) two different behaviours can be seen 
according to the two different regions of magnetism: the 
local moment magnetism (LMM) and the intermediate 
valence magnetism (IVM). For low e/, i.e. in the LMM 
regime (see inset of figure ^3|), the phase transition is of 
first order, for e/ in the IVM regime of second order. So 
the two distinct regions of ferromagnetism seen in figure 
^ show a different temperature dependent behaviour. 

In figure |l^, two interesting features can be ob- 
served. There can be a temperature driven increase in 
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Figure 12: s- and /-magnetization as function of the 
temperature T for ritot = 1-5, ?7 = 10 and F = 0.2. The 
temperature scale is explained in the text. 
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Figure 13: Same as figure |T^, but for ntot — 1-7. In the 
inset, the /-magnetization is shown for e/ = —0.3 (solid 
line) and e/ = —0.2 (dotted line). 



/-magnetization, and the magnetic polarization of the 
conduction band can change sign as function of temper- 
ature. 

A magnetization increasing with temperature seems 
to be quite extraordinary, and at first sight, unplausible 
for thermodynamic reasons. But on second sight, even a 
temperature-driven spontaneous ordering from a param- 
agnetic to a ferromagnetic state is possible In our 
case, the system does not undergo a phase transition, 
only the magnetization increases slightly for certain pa- 
rameters. This can be easily understood as a quasiparti- 
cle densities of states effect: If the system parameters are 
chosen, such that the lower edge of the upper hybridiza- 
tion spin |-band is close to, but still above the Fermi en- 
ergy, these states are unoccupied for T — 0. But already 
a small softening of the Fermi function can induce a 
strong increase in occupied spin-f states since this band 
edge generically consists of a rather sharp peak. The ra- 
tio of spin-t electrons to spin-| electrons increases, the 
magnetization rises. 

The other feature is observed in the lower picture 
of figure ^ showing the conduction band polarization. 
As already seen for T = 0, it can either be parallel 
(ferromagnetic coupling) or antiparallel (antiferromag- 
netic coupling) to the /-magnetization, but sometimes, 
e.g. for e/ = —0.05, it can also cross the iris = 0-axis 
at some temperature T*. Quite remarkable is the be- 
haviour of the /-magnetization which is continuous even 
at T = T*. This resembles the T = behaviour of the 
magnetization as function of ntot hi the IV regime dis- 
cussed in the previous section. Again, it is indicated that 
the coupling mechanism responsible for the alignment of 
the /-moments in the intermediate valence regime does 
not depend on the explicit polarization of the conduc- 
tion band. This polarization has to be understood as a 
consequence, and not the cause of the ordering of the 
/-system. 

In the figures |lj and [ij, the demagnetization be- 
haviour is examined by means of the quasi particle den- 
sities of states. These illustrate both the first order tran- 
sition in the LMM regime {cf — —0.2: figure |lj) and 
the second order transition in the IV region (e/ = —0.05: 
figure In both cases, the demagnetization is charac- 
terized by two effects, a bandshift and a weight transfer 
mechanism. The bandshift between spin | and spin J, 
band decreases with increasing temperature while at the 
same time a transfer of spectral weight between the hy- 
bridization bands occurs. For the spin |-bands, weight 
is transferred from the lower to the upper band, for the 
spin |-bands just the contrary takes place. 

For Ef = —0.2, ferromagnetism breaks down at 
T = Tc ~ 208, we see a first order phase transition. 
Contrarily, for e/ = —0.05, the demagntization behaves 
smoothly up to T = Tc « 114, thus a second order phase 
transition occurs. 

Finally, in figure the Curie-temperatures as func- 
tion of ntot (upper figure) and e / (lower figure) are pre- 
sented. The continuous vanishing of Tc as a function of 
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ritot corresponds to the second order phase transitions 
seen in figure |[ 

Comparing figure ^ and the Curie temperature as a 
function of ntot in the IV regime (figure |l6|) , one observes 
that the maximum oiT(.{ntot) coincides neither with the 
maximum of mf{ntot) not the maximum of mj at T = 0. 
Taking the Curie temperature as a measure of the ther- 
modynamic stabihty of ferromagnetism, this indicates 
again that the magnetic polarization of the conduction 
band does not directly cause the magnetic ordering of 
the /-moments. Thus, we conclude that the magnetic 
polarization of the conduction band is rather induced by 
the magnetically ordered atomic levels than responsible 
for the ordering itself. This supports our proposition 
that the magnetic ordering of the /-moments in the in- 
termediate valence regime is rather a single band effect 
within the /-system than caused by a spin-exchange be- 
tween /-levels and conduction band. 



Figure 14: /-electron quasiparticle density of state 
(QDOS) as a function of energy for different temper- 
atures (solid line: spin |, dashed line: spin |). The 
system parameters are ntot = 1-7, e/ = —0.2, U = 10 
and V = 0.2. The corresponding magnetization curve is 
presented in the inset of figure |l^. The system shows a 
first order phase transition at Tc « 208. 
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Figure 15: same as figure |T^, but for e-j — —0.05. This 
corresponds to the dashed line in figure 13, with Tc ~ 
114. 
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Figure 16: Curie-temperature Tc in the upper picture, 
as function of ntoti and in the lower as function of e/ 
{U = 10 and V = 0.2) 



In figure showing Tc as a function of ey , the split- 
ting of the two regions of ferromagnetism for ntot ~ 1-7 
can be seen again. The lowest boundary of ferromag- 
netism with respect to is given by the first order phase 
transition at the quantum critical point e|^^°j'^. Right 
above this critical point, both the Curie temperature as 
well as nif{T = 0) show their maximum. Here, in the 
LMM regime, the Schrieffer- Wolff transformation sug- 
gests that an effective coupling of /-moments is possible 
only via an (antiferro-) magnetic exchange with the con- 
duction band. But again, for the IVM region, no connec- 
tion can be made between the T = moments of the two 
electronic subsystems and the Curie temperature Tc- 
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4 Conclusions 

In this paper, we have presented a new approach to 
the periodic Anderson model. Our approximation is 
based on an exact mapping of the PAM onto an effec- 
tive Hubbard model. Since this effective model is clearly 
located in the strong coupling regime, a for this limit 
rather well motivated approximation, the SDA is used to 
find self-consistent solutions. Even though the particu- 
lar low-temperature properties of the Anderson model 
like a formation of a Kondo resonance are not accessible 
within our approximation, it is trustworthy concerning 
magnetism. A major advantage of this approach is its 
numerical simplicity, enabling calculations for the full 
parameter range and the accessibility of all major quan- 
tities of interest via the selfenergy given on the real axis. 

Using this method we have shown the T = phase 
diagram and have found two distinct regions of ferro- 
magnetism. The first of these is best described by the 
term "local moment magnetism" (LMM) while the other 
located within the intermediate valence region (IVM) 
displays itinerant character. 

The local moment region is close to the Kondo regime 
of the PAM. The magnetic moments within the /-system 
are fully polarized, the coupling between /-moments and 
conduction band is antiferromagnetic, as expected from 
the Schrieffer- Wolff transformation. The Curie temper- 
atures are rather high. 

On the other side, the magnetism in the interme- 
diate valence regime is quite weak, in a sense that the 
average magnetic moment per site is low and the Curie 
temperatures are smaller than in the LMM region. The 
magnetic polarization of the conduction band can cither 
be parallel or antiparallel to the /-moments, at some 
isolated points, the conduction band is not polarized at 
all. This means, the polarization of the conduction band 
should be of minor importance to the magnetic ordering 
of the /-moments. 

We conclude that the mechanism resposible for the 
ferromagnetic ordering is different for the LMM and 
the IVM regime. Whereas in the LMM regime, a spin- 
exchange according to the Schrieffer- Wolff transforma- 
tion causes ferromagnetic ordering of the /-moments, 
in the IV regime, the spin fluctuations still lead to a 
spin polarization of the conduction band, but rather 
more important for the ferromagnetic ordering in the 
/-levels are the charge fluctuations leading to an eff'ec- 
tive broadening of the /-Icivels to a band. This band 
has those properties expected to form ferromagnetism 
within a single-band Hubbard model. 
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